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Abstract. 12 CO and 13 CO J = 6-5 observations of 
IRC +10 216 with the JCMT are presented. A spheri- 
cally symmetric radiative transfer code is used to model 
these and other CO observations of the carbon star 
IRC +10 216/CW Leo. Compared to previously published 
model calculations a much larger set of observational data 
is used as constraints; on-source 12 CO and 13 CO J = 1-0 
up to J = 6-5 and mapping data taken with various tele- 
scopes, most of which are obtained from the literature. 

The gas temperature in the envelope is calculated tak- 
ing into account heating and cooling. The heating by dust- 
gas collisions and various other parameters (such as lumi- 
nosity over distance squared) are constrained from our 
previous modeling of the circumstellar dust shell. Photo- 
electric heating is taken into account. 

A grid of models is calculated with the following pa- 
rameters: luminosity (in the range 10 000 - 30 000, in steps 
of 5 000 L©), mass loss rate, dust-to-gas ratio, dust opacity 
and CO abundance. For each of the considered luminosi- 
ties a good fit to the on-source data can be obtained. A 
comparison with CO J = 1-0 data obtained from the lit- 
erature points towards a preferred luminosity of 10-15 000 

Notwithstanding the overall good agreement, there re- 
main discrepancies. The different observed 12 CO J = 3-2 
observations appear to be always larger than the model 
predictions. The observed 13 CO J = 6-5 is almost flat- 
topped, while the model gives a slight double-peaked pro- 
file. There is a large discrepancy with the single existing 
12 CO J = 7-6 observation. 

The best fitting models (for each of the considered lu- 
minosities) cannot accommodate the more extended emis- 
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sion seen in the mapping data. This is not due to an un- 
derestimate of the photoelectric effect. To fit the data for 
radial offsets >50" the mass loss rate must be a factor of 
5 higher in the outer envelope. Because the various sets of 
data for offsets >150" are not consistent with each other it 
is unclear if the enhancement in the mass loss rate extends 
beyond that radius. 

Visibility curves are calculated for comparison with fu- 
ture interferometric observations. These appear to be in- 
sensitive to luminosity and mass loss variations but should 
be good tracers of the geometry of the CO shell. 

A comparison is made between the mass loss rates 
and dust-to-gas ratios derived from the CO modelling and 
those derived from our previous dust modelling. To do 
this we make use of the relation M = rp (l — -p), 
which is valid for radiation pressure driven outflows. The 
best agreement is obtained for the model with 15 000 L . 
This agrees well with the luminosity range 7 700-12 500 
L/0 based on the period-luminosity relation for carbon mi- 
ras. 

In summary, we conclude that the likely luminosity of 
IRC +10 216 is between 10 000 and 15 000 L Q and that 
its distance is between 110 and 135 pc. The present-day 
mass loss rate is (1.5 ± 0.3) xl0~ 5 M Q yr _1 and the gas- 
to-dust ratio is a 700 ± 100. The dust opacity at 60 fim is 
found to be of order of 250 cm 2 gr _1 . The CO abundance 
is 1.1 x 10~ 3 relative to H 2 . 
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1. Introduction 

IRC +10 216 (= AFGL 1381 = CW Leo) is the best- 
studied carbon star, because it is intrinsically bright and 
nearby. One aspect of study has traditionally been the 
properties of the circumstellar shell: e.g., what is the mass 
loss rate, has the mass loss rate changed with time, what 
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kind of chemistry takes place, what is the geometry of the 
shell? 

In two recent papers we have made an extensive study 
of the circumstellar dust shell of IRC +10 216 (Groenewe- 
gen 1996, 1997, hereafter G97). A much larger set of obser- 
vational constraints was used than in any previous study. 
This allowed the derivation of the present-day dust mass 
loss rate, the ratio between the luminosity at maximum 
light divided by the distance squared, the dominant grain 
size, the effective temperature and the inner dust radius. 

In another paper we have mapped the 1.3 mm contin- 
uum emission around the star (Groenewegen et al. 1997). 
It was detected out to 50" at a noise level of 3a or bet- 
ter. There is evidence for phases of enhanced mass loss, 
but a quantitative analysis is hampered by the significant 
contribution of molecular lines to the continuum emission, 
and the uncertainty in the radial distribution of that con- 
tribution. 

Another way to study the circumstellar shell is by ob- 
serving the gas component. A good tracer of the (H2) gas 
is the CO molecule, and numerous CO observations ex- 
ist of IRC +10 216 (see the catalog of Loup et al. (1993) 
for a compilation) . A fair number of detailed models have 
been proposed to explain various CO observations. The 
pioneering work was done by Kwan & Hill (1977). For an 
assumed distance of 200 pc they deduced a mass loss rate 
of 2 x 10~ 5 M Q yr _1 , a CO abundance relative to H 2 of 
about 8 x 10~ 4 (both uncertain to a factor of 2), and a 
CO/ 13 CO isotope ratio of 35 ± 7. In a second paper Kwan 
& Linke (1982) deduced a mass loss rate of 4 x 10" 5 M 
yr^ 1 for a distance of 200 pc and a CO abundance of 6 x 
10~ 4 . Sahai (1987) finds, based on the modeling of high-J 
lines, that the temperature in the inner region is much 
larger than that predicted by the Kwan & Linkc model, 
and that the mass loss rate in the inner 6" region is 3.2 x 
10~ 5 M Q yr _1 and 50% higher in the extended envelope, 
assuming a distance of 300 pc and a CO abundance of 8 
x 10~ 4 . Truong-Bach et al. (1991) also find evidence for 
a hot inner core: a mass loss rate of 2.5 x 10~ 5 M Q yr _1 
in the inner 4.2" and a mass loss rate of 4 x 10~ 5 M 
yr _1 in the outer region (for a distance of 200 pc). Kast- 
ner (1992) found a best fit to some CO data for a distance 
of 150 pc and a mass loss rate of 2 x 10~ 5 M Q yr _1 . In 
a recent paper, Crosas & Menten (1997) also presented 
J = 6-5 12 CO and 13 CO data, and combined those with 
other data from the literature, and modelled the data. In 
this respect the philosophy behind that paper is similar 
to ours. The major differences with that paper are that 
the effective absorption coefficient is taken into account 
in a consistent way (see below and Sect. 7.2) and that the 
parameter space is scanned consistently no a-priori fixing 
any parameter. In Crosas & Menten, the dust-to-gas ratio 
is fixed at a value of 0.01. They find a good fit using a dis- 
tance of 150 pc, a mass loss rate of 3.25 x 10~ 5 M Q yr _1 , 
a CO abundance of 6 x 10~ 4 , and a dust-gas momentum 
transfer rate Q of 0.025 for the inner region (r < 10 16 cm), 



and 0.018 for the outer region. They conclude there is no 
need for a hot inner core. 

In all the model calculations mentioned above the ef- 
fective absorption coefficient Q is assumed to be a free 
parameter. Q describes the momentum transfer efficiency 
between the gas and the dust. However, we claim, and 
this is a key point in this paper, that effective absorption 
coefficient Q is not a free parameter but is related to the 
properties of the dust. This had been realised before but 
it requires the simultaneous fitting of the spectral energy 
distribution (SED) and the CO data in order to estimate 
Q in a self-consistent way. This approach was taken in the 
fitting of two OH/IR stars in Groenewegen (1994b) and is 
proposed here for IRC +10 216. 

The outline of the paper is as follows. In Sects. 2 and 
3 the model is briefly described. In Sect. 4 the 12 CO and 
13 CO J = 6-5 observations as well as the other CO data 
to which the models are compared are presented. In Sect. 
5 the results of the model calculations are described. So- 
called visibility curves, to be compared with future in- 
terferometric observations, are presented in Sect. 6. The 
paper ends with discussion and conclusions in Sect. 7. 



2. The model 

The calculations are performed with the molecular line 
emission radiative transfer model of Groenewegen (1994a; 
G94). This model is an improvement over and extension 
of the code by Morris et al. (1985) through the inclusion 
of a self-consistent determination of the gas temperature. 
The model was previously applied to two OH/IR stars at 
known distances (Groenewegen 1994b). 

The model is fully described in G94 and only some 
essential details are repeated here. The model assumes 
spherical symmetry and that the Sobolev approximation 
is valid to calculate the radiative transfer (this presumes 
that the local linewidth is much smaller than the expan- 
sion velocity). The line profiles are calculated without 
the Sobolev approximation using a ray-tracing method in- 
stead, assuming, in the present paper, a thermal velocity of 
0.5 km s _1 . The CO, respectively 13 CO or HCN, molecules 
are excited by: (1) collisions with H 2 and He molecules, 

(2) interaction with the 2.8 K background radiation, and 

(3) infrared radiation from a central blackbody of temper- 
ature Tbb and radius Rbb which leads to pumping from 
the v = into the v = 1 vibrational state. 

The heating rate by the photoelectric effect is discussed 
in G94. The heating rate per hydrogen molecule is not 
simply a constant as is often assumed, but the grain charge 
parameter is calculated for each radius and optical depth 
effects are taken into account. The photoelectric heating 
rate scales as (Go (Y/0.1)), with Go the UV flux of the 
diffuse interstellar medium, and Y the yield of dust grains. 
Go = 1 and Y = 0.1 arc used unless otherwise noted. 
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Molecule iransition 


Reference 


Telescope 


Beam width 
FWHM(") 




Remark 


CO 1-0 


Groenewegen et al. (1996) 


IRAM 


21.0 




primary constraint 




Truong-Bach et al. (1991) 


IRAM 


21.0 




map 




Groenewegen & Ludwig (1998) 


IRAM 


21.0 




map 




Knapp & Morris (1985) 


BTL 


100 








Huggins et al. (1988) 


NRAO 


50 




map 




Olofsson et al. (1993) 


SEST 


45 








Huggins et al. (1988) 


oso 


34 




map 


2-1 


Groenewegen et al. (1996) 


IRAM 


12.5 




primary constraint 




Truong-Bach et al. (1991) 


IRAM 


12.5 




map 




Groenewegen & Ludwig (1998) 


IRAM 


12.5 




map 




Huggins et al. (1988) 


NRAO 


31 




map 




Olofsson et al. (1993) 


SEST 


23 






3-2 


Groenewegen et al. (1996) 


JCMT 


14.3 




primary constraint 




Williams & White (1992) 


JCMT 


15 








Wang et al. (1994) 


CSO 


21.4 






4-3 


Williams & White (1992) 


JCMT 


11.0 




primary constraint 


6-5 


this paper 


JCMT 


55% in 7.0" 


+ 45% in 18" beam 


primary constraint 


7-6 


Wattenbach et al. (1988) 


UH 2.2m 


45 






CO 1-0 


Groenewegen et al. (1996) 


IRAM 


21.0 




primary constraint 




Knapp & Chang (1985) 


BTL 


105 






2-1 


Groenewegen et al. (1996) 


IRAM 


12.5 




primary constraint 


3-2 


Groenewegen et al. (1996) 


JCMT 


14.3 




primary constraint 




Williams & White (1992) 


JCMT 


15 








Wang et al. (1994) 


CSO 


20.3 






6-5 


this paper 


JCMT 


55% in 7.0" 


+ 45% in 18" beam 


primary constraint 



The heating rate per unit volume caused by dust-gas 
collisions is (see G94 for details): 



ap d \ M 

and the drift velocity in km s _1 given by: 
' L Q v{r 



Vdr 



1.4293 x 10" 



M 



0.5 



(1) 



(2) 



where pd is the dust grain density, a the grain size, L the 
stellar luminosity in solar units, M the mass loss rate in 
M Q /yr, Q the effective absorption coefficient (defined be- 
low in Eq. 4), v(j) the gas velocity in km s , fu e the 
number ratio of helium to hydrogen (a value of 0.1 is as- 
sumed) and the dust-to-gas ratio. 

The cooling rate is determined by molecular line cool- 
ing of 12 CO, 13 CO, HCN, H 2 and adiabatic cooling. The 
H 2 cooling rate is treated in an approximate way but is 
not very important (G94). The cooling rate due to 12 CO, 
13 CO and HCN is calculated from the molecular excita- 
tion code in an iterative process. For a given temperature 
structure the level populations are calculated, which al- 
lows the molecular cooling to be determined, and finally 
a new temperature structure is derived by solving the en- 
ergy balance equation (see G94). 



2.1. How does dust affect the molecular excitation model? 

As mentioned before dust grain-gas collisions is the most 
important heating mechanism and therefore the properties 
of the dust (size, opacity, dust-to-gas ratio) influence the 
gas temperature and hence the resulting CO line profiles. 

In G97 the spectral energy distribution (SED) and vis- 
ibility curves of IRC +10 216 were modelled using a dust 
radiative transfer (DRT) model and constraints were set 
to the dust properties. Table 2 lists for an assumed lumi- 
nosity, the derived mass loss rate and distance. 



Table 2. Default parameters 



Luminosity 


Distance 1 1 ' 


M dc{ 


(i) 


(L ) 


(PC) 


10" 5 


M Q yr- 1 


10 000 


110 


1.80 




15 000 


135 


2.20 




20 000 


156 


2.54 




25 000 


174 


2.84 




30 000 


191 


3.11 





^'Parameters that follow from the dust modeling in G97. The 



mass loss rate is based on <]/ = 0.005, kqo 
dust expansion velocity of 17.5 km s _1 . 



68 cm gr and a 
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Table 3. Grid of models ran 



* */1.2 


*/1.25 




tf/1.67 */1.8 */2 


*/3 


*/4 


*/5 


L = 10 000 L Q 














MdcfX 1.0 9 


- 


- 


10 


10 




9 


Ma c 1 1 9 








10 






(if, , / 1 A 








1 3 






M , , 1 1 Q 

JWdef / J-.o 






1 fi 








L = 15 000 L Q 














Af dcf x 1.2 


_ 


_ 


_ _ 


_ 


_ 


8 


M def x 1.0 9 


- 


- 


9 


9 


9 


- 


M, , 1 1 O in 
iw der / J-.z AU 






10 


]_]_ 






lif , , / U 19 

JV^der / J-.*± 1Z 






12 










14 












L = 20 000 L 














MdefX 1.2 


- 


- 






- 


8 


if, ,v 1 n q 
-^dcf «*> J-.u y 






1 fl 


Q 






iw def / 








12 






A^fj t / 1 4 






14 - - 








L = 25 000 L 














AfdefX 1.4 - 


- 


- 


7 


7 


7 


- 


MdefX 1.2 - 8 






8 




8 


8 


MdefX 1.0 9 


_ 


_ 


10 


9 






M dcf / 1.2 10 


- 


- 


10 


11 


- 


- 


Mj t / 1 4 19 

iw def / A.<± J-^ 






13 12 








!\> , , / 1 R 1 e; 
Jw def / r.U lO 


1 ^ 




12 








M, , / 1 s 17 














L = 30 000 L 














M de fX 1.2 - 8 














MdefX 1.0 10 






10 


10 






Afdef / 1.2 10 






11 








Mdef / 1.4 11 




12 











For every luminosity the combinations of models that have been calculated are indicated. The corresponding absolute dust 
opacity follows from Eq. (3). The default mass loss rate and the distance are different for each luminosity; see Table 2. The 
entry in the table is the best fitting CO abundance in units of 10~ 4 . Bold-faced numbers indicate best-fitting models discussed 
in Sect. 5. 



The shape of the calculated SED is determined com- 
pletely by the input spectrum of the underlying star, the 
inner dust radius (or equivalcntly, the temperature of the 
dust at the inner radius, T c ), and the optical depth. The 
optical depth is given by: 



ira 2 Q\ n d (r) dr 



M g Qx/a 
R* p A r c 



j; 



R(x) 
x 2 w{x) 



dx 



(3) 



where x — r/r c , M(r) — M R(x) and v(r) = Voo w(x). 
The normalized mass loss rate profile R(x) and the nor- 
malized velocity law w(x) should obey R(l) = 1 and w(oo) 
= 1, respectively. In the case of a constant mass loss rate 
and a constant velocity the integral in Eq. (3) is essen- 
tially unity since r max is always much larger than 1. The 
symbols and units in Eq. (3) are: the (present-day) mass 
loss rate M in M Q yr _1 , W the dust-to-gas mass ratio (as- 
sumed constant with radius), Q\ the extinction coefficient, 
a the grain size in cm (the model assumes a single grain 
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size), R± the stellar radius in solar radii, Voo the terminal 
velocity of the dust in km s _1 , pa the dust grain specific 



density in g cm 



the inner dust radius in units of 



stellar radii and r max the outer radius in units of r c . The 
mass loss rates derived in G97 are based on an assumed 
dust velocity of 17.5 km s _1 , an absolute opacity at 60 fim 
of 68 cm 2 g _1 and a dust-to-gas ratio of 0.005. 

The relation between the wavelength dependent ab- 
sorption coefficient Q\ in Eq. (3) and the effective ab- 
sorption coefficient Q in Eq. (1) is: 



Q = 



JF X Q X d\ 
J'F x dX 



(4) 



where F\ is the emerging flux from the central star and the 
dust shell at infinity. In principle Q depends on the radius 
since F\ is continuously changed by the dust emission but 
this is only important near the inner radius where the 
dust temperature changes rapidly (see Habing et al. 1994). 
The molecular transitions under discussion here originate 
from several tens to hundreds of stellar radii where F\ is 
essentially independent of the radial distance. 

Radiative pumping of molecules is provided by ther- 
mal emission from hot dust close to the star. In molecular 
emission models this is usually represented by a blackbody 
of temperature Tbb and radius -Rbb • These quantities can 
be estimated from the DRT-models. At each gridpoint in 
the DRT-model the blackbody temperature of the radia- 
tion field is determined. In this way a realistic estimate of 
Tbb and i?BB is obtained. 

In G97 we showed that the near-infrared visibility 
curves are very sensitive to the grain size and that a main 
grain size of 0.16 ± 0.01 pm could fit the visibility curves. 
Therefore a grain size of 0.16 pm is adopted in all calcu- 
lations. 

Finally, DRT-models provide To.i, the optical depth 
at 1000 A which is needed to calculate the photoelectric 
heating (see G94 for details) and the dust temperature 
profile which is needed in the heating rate due to the gas- 
dust temperature difference (see G94). 

2.2. The outer radius and photodissociation 

The outer radius of the CO shell is determined by pho- 
todissociation. Mamon et al. (1988) have investigated this 
effect and found that the CO abundance as function of 
radius, relative to the value close to the star, can be ap- 
proximated as: 



Xco 



-ln(2) (r/r 1/2 )° 



(5) 



where r 1 / 2 and a depend on M and v and are tabulated by 
Mamon et al. (1988). Here we instead use the analytical 
fit to these results by Stanek et al. (1995): 



+ 7.5 x 10 15 ^— J cm 
and Kwan & Webster (1993) 

0.09 

/ M 1 5 \ 

a = 2.79 



M 15 
10" 5 ~ 



(6) 



(7) 



where fco is the number ratio of CO to H2 molecules 
close to the star, M is the mass loss rate in M Q yr _1 and 
v the gas velocity in km s _1 . The photodissociation radius 
of 13 CO is to first approximation assumed to be equal to 
that of 12 CO (see Mamon et al. 1988). The outer radius of 
the CO shell is set at the radius where the CO abundance 
has dropped to 1% in the calculations. Equations. (6) and 
(7) are used with a gas velocity of 14.5 km s _1 , which is 
the observed terminal velocity of the gas. 

3. Best guesses for the dust-to-gas ratio and CO 
abundance 

Two important free parameters are the dust-to-gas ratio 
and the abundance of CO molecules relative to H2. Both 
quantities can not take arbitrary values. 

There is a physical upper limit to the dust-to-gas ratio 
based on the number of atoms that can condense into 
dust. Using the continuity equation for the gas and the 
dust and assuming that the dust is 100% carbonaceous 
one finds that the theoretical dust-to-gas ratio is given 
by: 



tt = / c (C/Q-l)gg — " sas + " dr 

n H 1.4 v s 



(8) 



'''gas 



n/2 



5.4 x 10 16 



M 

lO- 6 



0.65 



/. 



CO 



-0.55 

15 J { 8 x 10- 4 



0.55 



where f c is the degree of condensation of the dust and 
C/O is the number ratio of carbon to oxygen atoms in 
the gas phase. Taking the observed gas velocity of 14.5 
km s _1 , and assuming a drift velocity of 3.0 km s -1 and 
a cosmic oxygen abundance of 8.70 (on a scale where H = 
12.0) then the theoretically predicted dust-to-gas ratio is 
# = 5.18 x 10- 3 / c (C/O -1). 

For disk planetary nebulae it found that C/O < 4 (see 
e.g. Groenewegen & de Jong 1994 and references therein) 
and because / c < 1 by definition, it follows that ^ must be 
smaller than 0.016. This upper limit is valid for all carbon 
stars. Calculations show that f c is about 0.4 (within a fac- 
tor of 2) (see Fleischer et al. 1995, Winters et al. 1994a,b) 
and for IRC +10 216 a C/O ratio of 2-3 seems appropri- 
ate (G97). In this case * is in the range 0.001 - 0.008. In 
the model calculations we use dust-to-gas ratios between 
0.001 and 0.005, in agreement with the theoretically ex- 
pected range. 

For fco, defined as the n(CO)/n(H 2 ) ratio close to 
the star, one may assume that in a carbon star all oxygen 
atoms are locked up in CO. Using an oxygen abundance 
of 8.70 ± 0.1 it immediately follows that f co = (10 ± 2) x 
10~ 4 . All the best fitting models have derived CO abun- 
dances in this range. Note that some of the other models 
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Table 4. Fit to secondary constraints 



Transition Telescope Observed Peak Model Peak Temperature(K) 







Temp. (K) 


L = 10 000 L 


15 000 


20 000 


25 000 


30 000 


i2 CO (1-0) 


BTL 


5.1 


4.9 


5.5 


5.9 


5.5 


5.8 




NRAO 


9.6 


9.3 


10.3 


10.9 


10.3 


10.8 




SEST 


10.7 


10.2 


11.2 


11.9 


11.2 


11.7 




oso 


10.0 


13.2 


14.7 


14.9 


14.2 


14.5 


12 CO(2-l) 


NRAO 


21.5 


15.7 


16.8 


17.5 


16.6 


17.2 




SEST 


25.4 


21.1 


22.1 


23.0 


21.9 


22.7 


12 CO(3-2) 


CSO 


32.5 




22.6 








12 CO(7-6) 


UH2.2 


9.0 


2.87 


2.66 


2.78 


2.65 


2.75 


13 CO(1-0) 


BTL 


0.38 


0.27 


0.29 


0.32 


0.30 


0.32 


13 CO(3-2) 


CSO 


2.9(0) 




2.62 









Wang et al. (1994) give a value ol 3.9 K, but this apparently refers to the peak temperature of the horn, in this double-peaked 
profile. The observed temperature at the line center is read of their Fig. 2 and then appropriately scaled to give the main-beam 
temperature. 



discussed in the literature (Kwan & Linke 1982, Crosas & 
Menten 1997) find a CO abundance of 6 x 1CT 4 , which is 
in disagreement with the theoretically expected value. 

4. The CO data 

A wealth of CO data exist for IRC +10 216. In Table 1 we 
list the data that are used as constraints here. Except for 
the 12 CO and 13 CO J = 6-5 spectra, all data are taken 
from the literature. 

The 12 CO and 13 CO J = 6-5 observations were carried 
out using the James Clerk Maxwell Telescope (JCMT), lo- 
cated on Mauna Kea, Hawaii. The detector was an SIS re- 
ceiver (RxG) available to JCMT via a collaborative agree- 
ment with the Max-Planck-Institut fur Extraterrestischc 
Physik in Garching (Germany). A gunn oscillator yields 
continuous coverage across the 650-692 GHz band, which 
includes the 12 CO and 13 CO J = 6-5 lines. The backend 
was an acousto-optic spectrometer. Before the observa- 
tions were made, the planets Jupiter and Mars were used 
to determine the beamshape during good weather with a 
30% zenith transmission. It was found that, although the 
JCMT beam is complex at these high frequencies, it can be 
well described by a two-component beam with 55% of the 
power in the 7" FWHM diffraction component and 45% 
of the power in a 18" FWHM beam. The overall beam ef- 
ficiency in this compact beam is 27% and was determined 
by taking spectra of Jupiter and Mars. When taking all 
errors into account we expect the calibration error to be 
less than 20%. 

The 12 CO J = 6-5 observations were made on April 21 
and 25, 1995. The pointing was done using the continuum 
receiver UKT14. The relative pointing difference between 
the RxG and UKT14 receivers was determined using very 
bright sources and was accurately known. After the tele- 
scope was successfully pointed at the source using UKT14, 
a switch was made to RxG and the relative pointing cor- 



rection was applied. During the observations IRC +10 216 
was always above 70 degrees elevation. The zenith sky 
transmission was 17% on April 21 and 19% on April 25 
yielding system temperatures of 17 000 and 22 000 K re- 
spectively. Four scans of 4 min each were made, three of 
which were on April 21, for a total integration time of 16 
min. On April 21 the on-source spectra were part of a 5- 
point map with a beam separation of 8"; on April 25 the 
on-source spectrum was part of a 5-point map with a 10" 
beam separation. 

In Figs. 1-5 one can notice a bump in the spectrum 
centered near a velocity of km s _1 . It is not present 
in the J = 6-5 spectrum of Crosas & Menten (1997). We 
checked the JPL line catalog and no suitable candidate 
could be found. We also checked for lines coming in from 
the sideband with the same result. In any case this uniden- 
tified line does not influence any conclusions regarding the 
(fitting of the) J = 6-5 profile. 

The 13 CO J = 6-5 observations were made on April 30, 
1995. The pointing was done as described above for the 
12 CO observations. During the observations IRC +10 216 
was between 65 and 70 degrees elevation. The zenith sky 
transmission was 32% yielding a system temperature of 
10 000 K. Seven scans of 8 min each were made for a total 
integration time of 56 min. No mapping was done in this 
line. 

The above data were added after visual inspection of 
the individual scans and base line subtraction. The rms 
noise in the final 12 CO J = 6-5 spectrum is 2.4 K while 
the line is about 45 K. The rms noise in the final 13 CO J 
= 6-5 spectrum is 0.8 K while the line is about 7 K. 

The primary constraints in this paper are the J = 6-5 
profiles discussed above, the 12 CO J = 4-3 profiles from 
Williams & White (1992) and 12 CO and 13 CO J = 1-0, 
2-1 and 3-2 profiles observed at IRAM and JCMT and 
discussed in Groenewegen et al. (1996). In addition a se- 
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lection of other data was taken from the literature in par- 
ticular that taken with smaller telescopes. 
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Fig. 1. Model with L — 10 000 L©, and default parameters for 
M, and Q. Upper left panel indicates gas temperature (solid 
line) and excitation temperature of the 12 CO (1-0) transition 
(dashed line). The rise in the gas temperature at large distances 
from the star is due to photoelectric heating. Other panels 
indicate 12 CO (upper part) and 13 CO (lower four panels) line 
profiles. In each case the histogram indicates the observations. 
They are the "primary constraints" of Table 1. The solid line 
indicates the model, while the dashed line represents the model 
scaled to the observed maximum, in order to better compare 
the line shapes. 



5. Model calculations 

A large grid of models was calculated, varying the lumi- 
nosity (and the distance correspondingly), the mass loss 
rate, the dust-to-gas ratio, the absolute dust opacity and 
fco while fixing fco/f 13 co a t a value of 44, a value ac- 
curately known from optically thin molecular lines (see 
Kahanc ct al. 1992). The velocity law of the gas is param- 
eterised by v(r) = 14.5 km s" 1 (1 - 2.73 x 10 14 cm (d/135 
pc)/r) - 5 . 

From the dust modeling in G97 the following parame- 
ters were derived: Q = 0.0419; a = 0.16 /zm, inner (dust) 
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Fig. 2. Model with L = 30 000 L©, and default parameters 
for M, * and Q. 

radius = 3.17 x 10 14 cm (d/135 pc); T BB = 860 K; R BB = 
6.18 x 10 14 cm (d/135 pc); r .i = 1.84 x 10 16 cm (d/135 
pc)/r, with r the radial distance to the star. Given the 
observed integrated flux, Table 2 lists the distance corre- 
sponding to several chosen values of the luminosity, and 
the mass loss rate following from the dust model. 

The actual combinations of parameters used in the var- 
ious CO model calculations are listed in Table 3. For ex- 
ample, take the entry with L = 20 000 L Q , "M dof /1.4" and 
"\Jj/1.67" . The default mass loss rate for that luminosity is 
2.54 x 10 -5 M yr -1 , and the default dust-to-gas ratio is 
0.005 (see Table 2). This therefore means that the model 
was run with M= 1.81 x 10~ 5 M Q yr" 1 and * = 0.0030. 
From Eq. (3) it follows that the dust opacity was (1.4 x 
1.67) x 68 = 160 cm 2 gr _1 . For every model combination 
the CO abundance ratio was tuned to give the best pos- 
sible fit to both the 12 CO and 13 CO line profiles marked 
as "primary constraint" in Table 1. The 13 CO abundance 
was fixed at (/co/44) as explained above. 

Figures 1-2 present the fits to the "primary con- 
straints" for models with the default values for the mass 
loss rate, opacity and dust-to-gas ratio for the two extreme 
luminosities. Larger luminosities result in higher predicted 
brightness temperatures. This is most clearly seen in the 
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Fig. 3. Best fitting model with L = 10 000 L . 



Fig. 4. Best fitting model with L — 15 000 L . 



12 CO(1-0) profile. Overall however, the effect of luminos- 
ity appears not be very significant, and therefore only the 
extreme models with L = 10 000 and 30 000 are shown. It 
is also clear that these models (and the models with inter- 
mediate luminosities) fail to fit the "primary constraints" . 

Figures 3-5 present the best fitting models for three 
different luminosities. These models have a bold-faced en- 
try in Table 3. For every luminosity a model can be found 
that fits the "primary constraints" almost equally well. 
The high-luminosity models arc not shown because they 
fit qualitatively equally good, and such high luminosities 
can be ruled out on other grounds (see later). In all cases 
the models predict brightness temperature for the 12 CO 
(3-2) line that is too low. This may indicate a systematic 
uncertainty in the calibration (also see discussion in Sect. 
5.1.1). Also, all the models predict a 12 CO J = 1-0 line 
that is too bright. This difference is lowest for the 10 000 
Lq model, where the difference is well within the calibra- 
tion uncertainty. The remaining discrepancy is the 13 CO 
(6-5) line. The absolute intensity is well reproduced but 
the line shape does not fit very well. None of the mod- 
els calculated predict a flat-topped 13 CO (6-5) line. Also 
the recent model calculation of Crosas & Menten (1997) 



predicts a slightly double-peaked profile. Note that their 
model does not reproduce the absolute intensity neither. 

It was verified that this discrepancy is not due to a 
pointing off-set. A model calculated 3.5" off-centre still 
produced a slight double-peaked profile. A model at 7" 
off-set produced at flat-topped profile but the absolute 
intensity was a factor of 3 below the observed one. 

This discrepancy hints to the fact that possibly one 
of the fundamental assumptions in the model may no 
longer be full-filled: spherical symmetry or homogeneity. 
This can in principle be addressed by high spatial resolu- 
tion, high frequency mapping observations. Problems as- 
sociated with these high frequency observations are that 
the beam is usually small, and random pointing errors can 
already be a significant fraction of the primary beam, and 
the influence of the error beam, which makes interpreta- 
tion not straightforward. 

5.1. The other constraints 
5.1.1. On-source data 

We compared the best fitting model at each luminosity 
with the line intensities and line profiles for all on-source 
observations that were not used as primary constraints 



M.A.T. Groenewegen et al.: IRC +10 216 revisited II: the circumstellar CO shell 



9 




12 3 -40 -20 

log radius [10 15 cm] V LgR [km/s] 




-40 -20 -40 -20 




-40 -20 -40 -20 



LP 


• 


1-0 ; 





-40 -20 -40 -20 




-40 -20 -40 -20 

V IHi [km/s] V LSR [km/s] 



Fig. 5. Best fitting model with L = 20 000 L Q . 

(see Table 1). The peak intensities predicted by the model 
are listed in Table 4, together with the observational data 
from the respective papers. 

The reasoning behind this exercise is that these other 
J = TO and 2-1 observations were performed with smaller 
telescopes with large beam sizes, and so these observations 
sample a greater extent of the envelope than the higher 
resolution IRAM observations. 

The J = 1-0 BTL, NRAO and SEST data are equally 
well fitted by the five models, given the typical la uncer- 
tainties of 10% in the observations. In fact, all three obser- 
vations would point to a luminosity between 10 000 and 
15 000 L Q . Regarding the OSO observations there must 
have been a calibration error. Indeed, Huggins et al. men- 
tion that no beam efficiency correction is applied and that 
the observations have been scaled to an earlier 1982 OSO 
observation. A more recent OSO observation (Olofsson et 
al. 1993) indicates a peak temperature of 12.5 K which is 
low but formally in agreement with the models. The 13 CO 
J = 1-0 BTL observation agrees reasonably well with the 
models. 

The J = 2-1 NRAO observations are barely in agree- 
ment with the models. Other NRAO observations find 
similarly high on-source temperatures (Wannier & Sahai 



1986, Wannier et al. 1990). The SEST data agrees with 
the models. 

With respect to the J = 7-6 observations there is a 
discrepancy of more than a factor of three. We have no 
explanation for this. This is the only existing 7-6 obser- 
vation of this star, and we suggest that this transition be 
re-observed. As the J = 6-5 data is well fitted, there is no 
reason to believe that our model is in error by such a large 
factor. 

The referee kindly pointed us to the paper of Wang 
et al. (1994) which was not considered originally and dis- 
cusses 12 CO and 13 CO (3-2) data taken with the CSO. For 
one model, the relevant temperatures were calculated, and 
are included in Table 4. The 13 CO (3-2) model prediction 
is in excellent agreement with the observation, while the 
12 CO (3-2) model prediction is significantly lower. 

Williams & White (1992) also present 12 CO and 13 CO 
(3-2) data taken with the JCMT, and can therefore be 
directly compared to our observations. For the 12 CO line 
they find a main-beam temperature of 32 K at 2" off- 
set, and estimated an on-source temperature of about 37 
K. We find 41 K, in excellent agreement. For 13 CO they 
find a line-center temperature of 3 K (at a position ±5"), 
while we found 3.8 K. This agreement is acceptable, and 
it should be pointed out that their spectrum has a consid- 
erable lower S/N than ours. 

In conclusion, the two sets of JCMT data are in agree- 
ment with each other. For both CSO and JCMT the pre- 
dicted 12 CO (3-2) line is less intense than observed. The 
same is found by Crosas & Menten (1997) who compare 
their model to the Wang et al. data and find a model tem- 
perature in agreement with our value in Table 4 and hence 
lower than observed. In Sect. 5.0 it was suggested that this 
might be due to a calibration problem. This would imply 
calibration problems at both CSO and JCMT, which ap- 
pears unlikely. On the other hand a physical cause also 
appears unlikely. As temperature and density are smooth 
functions of radius, the profiles are also expected to be 
smooth function of increasing J-level. Since the 12 CO J=2- 
1 and 4-3 levels are well fitted, there is no a-priori reason 
to believe that the 3-2 level would be off. One could as- 
sume that at a radius, where the emission from the 3-2 
level is strongest, the density is lower than predicted from 
the standard model, but this would also then imply a dif- 
ferent 12 CO/ 13 CO ratio at that particular radius only as 
the the 13 CO 3-2 profile is predicted correctly. 

5.1.2. Mapping data 

In Fig. 6 the four best fitting models are compared to the 
mapping data of Huggins et al. (1988), Truong-Bach et 
al. (1991) and Groenewegen & Ludwig (1998), both on a 
linear and logarithmic scale. For the IRAM observations 
integrated intensities are compared, while for the OSO and 
NRAO observations by Huggins et al. only peak intensities 
are given in their paper. There exists at map in the 12 CO 
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Fig. 6. Comparison between observations (symbols with error bars) and models (lines) as a function of radial offset. IRAM 
observations come from Truong-Bach et al. (1991; squares), and Groenewegen & Ludwig (1998; triangles). The OSO and NRAO 
observations are from Huggins et al. (1988). Best fitting models are for L = 10 000 (solid), 15 000 (dashed), 25 000 (dotted), 
30 000 Lq (dash-dot-dot-dotted), and default values for photoelectric heating and constant mass loss rate. For the IRAM 
observations normalised integrated intensities, and for the OSO/NRAO observations normalised peak brightness temperatures 
are shown, on both a linear (left hand) and logarithmic scale (right hand). 
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J = 3-2 line (Williams & White 1992), but they find that 
the emission appears not to be extended with respect to 
the beam out to their maximum off-set of 50". Therefore 
we did not consider this data. 

A first remark is that there is extremely good agree- 
ment between the two sets of IRAM observations. The 
error bars plotted indicate the variation in the ratio of the 
integrated intensity relative to the on-source value at the 
different map positions for a given radial offset. 

At first glance there appears to be good agreement 
between the models and observations when plotted on a 
linear scale, but, as was noted by Huggins et al. (1988), 
clear discrepancies are seen when plotted on a logarith- 
mic scale. Although the absolute calibration uncertainties 
are of order 10% (and larger in the case of the OSO J = 
1-0 observation), the relative intensities are much better 
determined and only limited by noise. 

The models and observations start to deviate at ap- 
proximately 50" offset from the center position, the exact 
value depending on the beam-size. There are two possible 
explanations for this. First, the mass loss rate may have 
been higher in the past. Second, the photoelectric effect 
may be more important than assumed. 

Models with an increased mass loss in the past have 
been constructed to fit the data for luminosities of 15 000 
and 30 000 L . A good fit is obtained when the mass loss 
rate is higher by a factor of 5 for radial offsets > 50". These 
models are shown in Fig. 7 on a logarithmic scale only. The 
fit to the IRAM data is very good. That to the NRAO J 
= 2-1 and OSO J = 1-0 data is acceptable. On the other 
hand there is an obvious discrepancy between the OSO 
and NRAO J = 1-0 data. In the present model we as- 
sumed that the increase in the mass loss rate is constant 
and extends to the outer radius. Given the presently avail- 
able relatively poor data for radial offsets >90", and the 
discrepancies for radial offsets >150", it can not be ex- 
cluded that the mass loss rate is lower again in the outer 
parts of the envelope. The error in the radius of the on-set 
of the enhanced mass loss rate and the factor by which 
this is so are approximately 10" and less than a factor of 
2, respectively. From the radial dependence observed with 
IRAM compared to the constant mass loss rate case it is 
immediately clear that there is a break in the observa- 
tions between 40 and 60". It was verified numerically that 
a model with an increase in mass loss rate of a factor of 
3 did not fit the data. Although we do not claim that our 
best-fitting model is unique, this gives some constraints 
on the associated uncertainties of this model. 

As the increase in the mass loss rate occurs at off- 
sets much larger than the respective beam sizes, negligi- 
ble changes are expected in the on-source fluxes. This was 
verified numerically. 

Models with increased photoelectric heating were also 
considered (for the case of a constant mass loss rate and L 
= 15 000 L©). Even a model with the parameter (G x Y; 
see Sect. 2) a factor six larger than the default value could 



not fit the IRAM data, although it could fit some of the 
NRAO/OSO data. The reason is, that for the photoelec- 
tric effect to be important even at 50" from the star the 
interstellar radiation field must be very strong. Apart from 
the fact that the model with increased photoelectric effect 
does not fit the IRAM data, such a model appears also 
unphysical, as for a distance of 150 pc and its galactic lat- 
itude IRC +10 216 is at 110 pc from the galactic plane, 
where one does not expect an increase in the diffuse inter- 
stellar radiation field by such a large amount. 

6. Visibility curves 

In principle, intcrfcrometric observations yield additional 
information on the structure of the circumstellar envelope. 
In an article in the JCMT newsletter, Hills (1996) briefly 
mentions that interferometric observations at the wave- 
length of the J = 4-3 line had been performed with the 
JCMT-CSO interferometer. To provide comparison with 
future observations, we have calculated the so-called visi- 
bility curves at the center wavelengths of the J = 1-0 up 
to 6-5 12 CO transitions for the best fitting models with 
constant mass loss rate. For a brief summary on what a 
visibility curve is, and how they are evaluated, we refer to 
Sect. 2.1 of G97. The results of the calculations are plotted 
in Fig. 8. 

The model with mass loss variations as described above 
yield results indistinguishable from the constant mass loss 
rate case. This is due to the fact that these mass loss 
variations occur at large radii. 

At the same time this implies that the visibility curves 
may be most sensitive to the geometry of the CO shell. 
Asymmetries would be revealed by obtaining the visibility 
curves for different position angles. 

7. Discussion 

7.1. General comments 

We have made the most consistent modelling of the CO 
and dust envelope around IRC +10 216 so far. In par- 
ticular we did not assume that the gas-dust momentum 
transfer rate is a free parameter but constrained it from 
our previous dust modelling. The physical range of the 
dust-to-gas ratio and the CO abundance are considered, 
contrary to some other models that find unrealistically low 
values for the CO abundance. 

Contrary to Sahai (1987) and Truong-Bach et al. 
(1991), and in agreement with Crosas & Menten (1997), 
we find no evidence for a hot core in the inner 6" region. 
Specifically with respect to the Truong-Bach et al. result, 
Crosas & Menten suggest that they took too few rotational 
levels into account and therefore underestimated the rota- 
tional cooling in the inner part. Truong-Bach et al. infer- 
ence was based on an old J = 6-5 observation (Koepf et 
al. 1982). For the best-fitting 15 000 L model we find a 
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Fig. 8. Visibility curves for the best fitting models with constant mass loss rate, L = 10 000 (solid), 15 000 (dashed), 20 000 
(dot-dashed), 25 000 (dotted), 30 000 L (dash-dot-dot-dotted), for different transitions. The calculations are performed with 
a resolution of 0.01 arcsec -1 . 



peak temperature of 5.6 K in the corresponding 35" beam, 
almost a factor of three below the reported value. 

Cernicharo et al. (1996) discuss an ISO LWS grating 
observation of IRC +10 216. On top of the dust contin- 
uum there is a forest of HCN and CO lines. They infer 
that temperatures of order 700-1500 K are needed to ex- 
plain these observations. Specifically they adopt a temper- 
ature of 1200 K for radii smaller than 6 x 10 14 cm. From 
our self-consistent temperature determination we find a 
gas temperature at the inner radius (3.2 x 10 cm for a 
distance of 135 pc) of 1100 K. This agreement is encour- 
aging. 

In Groenewegen et al. (1997) we mapped IRC +10 216 
in the 1.3 mm continuum and found evidence for phases 
of enhanced mass loss in the past at radial distances of 
5 and 20" from the star. Although we find evidence for 
a higher mass loss rate in the past from the CO observa- 
tions, the radial scales are different. This is not necessarily 
incompatible. First, the largest variation found from the 



continuum map is only a factor of two at most and takes 
place from 10 to 30" at most. This may be not detectable 
in the CO. Secondly, as the continuum emission is detected 
at a significant level only out to 45" it is not possibly to 
prove or disprove if the variations seen in CO are also seen 
in the continuum emission. 

7.2. The influence of the radial dependence of Q 

Crosas & Menten (1997) find that they require a higher 
value for Q in the inner region, specifically Q = 0.025 
for r <,10 16 cm (that is 20 times their inner radius) and 
Q = 0.018 for r > 10 16 cm. This touches upon the ex- 
pectation (see Sect. 2.1) that Q should be a function of 
distance to the star because the effective wavelength of 
the emerging spectrum is shifting to longer wavelength 
due to the influence of dust. Crosas & Menten claim that 
their adopted behaviour for Q is in agreement with recent 
calculations, e.g. Habing et al. (1994), which showed that 
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Fig. 7. As previous figure but now for models including a 
higher mass loss rate in the past, as discussed in the text. L = 
15 000 (dashed), 30 000 L (dash-dot-dot-dotted). 
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Fig. 9. Radial dependence of the gas-dust momentum transfer 
efficiency Q. 

Q(r) is a strongly decreasing function of r. We in fact 
took the results by Habing et al. to assume a constant Q 
throughout the envelope. We will now address the influ- 
ence of this assumption by taking the radial dependence of 
Q into account. Figure 9 shows the radial dependence of 
Q as calculated from the best-fitting dust model. It shows 
that Q is rapidly dropping from large values to a value of 
0.0419 adopted in all calculation so far. However, and con- 
trary to the argument used by Crosas & Menten, it is clear 
that Q starts to deviate from its value at infinity only for 
radial distances much smaller than 10 16 cm. In Table 5 on- 
source peak temperatures for various 12 CO transitions are 
collected for the best- fitting model with L = 15 000 L Q , 
and different assumptions about Q. The list includes the 
J = 9-8 and 12-11 transitions which can not be measured 
from the ground, but can be from space. Shown are the 
standard case with a constant Q down to the inner radius, 
and various models where the radial dependence of Q is 
taken into account down to a certain inner radius. For in- 
ner radii smaller than 5 x 10 14 cm the model experiences 
some convergence problems. Since a greater value for Q 
means more heating (see Eq. 1) it is expected that models 
which take into account the radial dependence of Q lead 
to higher temperatures in the inner region and hence more 
emission from the high-J transitions. 

The result of the calculations show two things. For the 
low- J transitions, where the radial dependence of Q should 
be of no importance since the excitation temperatures are 
low, it is important to choose the inner radius not too big 
as this will result in an underestimate of the brightness 
temperature. This might also be part of the explanation 
for the need of an inner hot core by Truong-Bach et al. 
They assumed an inner radius of 2" or 6.0 x 10 15 cm for 
their adopted distance. This results in an underestimate 
of the true brightness temperature. If they had chosen a 
smaller inner radius (equal to the inner dust radius, which 
is in fact a physically plausible choice) the need, if any, for 
a hotter inner core would have been less. 

For the J = 3-2 and 4-3 transition the effects are 
smaller than 5% in the present case. The J = 7-6 result 
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Table 5. Influence of the radial dependence of Q. On-source peak temperatures for various transitions. 



Transition 


Beam 


r inn = 3.2 xlO 14 cm 


3.2 xl0 lb 


1.7 xl0 ib 


7.9 xlO 14 


5.0 xlO 14 




(") 


Q(n„„) = 0.0419 


0.0466 


0.0503 


0.0629 


0.0838 
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21 


T mb = 19.2 
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37.4 


32.9 


36.2 
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38.2 


3-2 


14.3 
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29.4 


32.1 


34.0 


34.4 


4-3 
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41.4 


34.5 
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42.5 


43.6 


6-5 
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62.0 


45.0 
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63.8 


67.3 


7-6 


45 


2.7 


2.4 


2.6 


2.8 


2.9 


9-8 


80 


0.53 


0.43 


0.49 


0.56 


0.59 


12-11 


60 


0.55 


0.37 


0.46 


0.58 


0.63 



shows that the effect becomes important for higher levels, 
and smaller beam sizes. As the high-J levels are generally 
measured with large beams, the effects are again moderate 
(of order 10%). 

7.3. The overall best fit 

Can a best model be identified? Not from the dust mod- 
elling alone as is mentioned in G97, since the shape of the 
spectral energy distribution is determined by the optical 
depth and the luminosity is a scaling parameter (see Ta- 
ble 2). The reason to consider CO observations as well, is 
that the heating rate is a non-linear function of luminos- 
ity (see Eq. 1). It turns out that for a plausible range of 
luminosities an almost equally best-fitting model to the 
primary constraints can be found. 

Some of the published CO J = 1-0 observations point 
to a luminosity <,15 000 L©. An independent constraint 
is the period-luminosity relation for carbon miras (Groe- 
newegen & Whitelock 1996). For a period of 649 days (Le 
Bertre 1992) it implies a luminosity of 9850 L©. The un- 
certainty in the period of about 10 days implies an uncer- 
tainty of about 150 L©. The uncertainty in the fit itself is 
0.26 in mboi which implies an error of 30% in luminosity. 
Thus from the PL-relation one derives therefore a lumi- 
nosity between 7 700 and 12 500 L©. 

An additional constraint comes from the observed ter- 
minal gas velocity. The assumption that radiation pres- 
sure on dust drives the outflow puts constraints on the 
mass loss rate and other parameters. Basically the formula 
M (vqo - Do) = t f f (l — must be obeyed (see Netzer 
& Elitzur 1993, Ivezic & Elitzur 1995, Groenewegen et al. 
1998), with v the gas velocity at the inner radius, tf the 
flux-weighted optical depth (known from the DRT mod- 
elling), and r is the ratio of radiation pressure to gravita- 
tional pull which scales like V ~ Q - ?. ^ — . Details 

1 a p d M„ "Uoo+Ddrift 

on how this system of equations is solved iteratively is ex- 
plained in Groenewegen et al. (1998). The equation allows 
one to verify that the mass loss rate and dust-to-gas ratio 
derived from the CO modelling are consistent with those 
derived from the DRT modelling using the above equation 



and procedure. Unknowns are the mass of the star (M* = 
0.8 M adopted), the gas velocity at the inner radius (u 
= 1 km s _1 adopted). The dust opacity is the one given by 
the best-fitting CO model for each luminosity. The results 
are listed in Table 6. 



Table 6. Comparison between CO and infrared modelling 



Lum. 




Mco 






Mir 




*IR 


(Le) 


cm /gr 


M Q /yr 




10~ 2 


Mo/yr 




10~ 2 


10 000 


326 


1.50x 10" 




0.125 


l.lOx 10" 


b 


0.137 


15 000 


245 


1.83x 10" 


-5 


0.167 


1.73x 10" 


-5 


0.135 


20 000 


245 


2.11x 10" 


-5 


0.167 


2.37x 10" 


-5 


0.113 


25 000 


283 


3.40x 10" 


-5 


0.100 


3.02x 10" 


-5 


0.099 


30 000 


204 


3.11x 10" 


-5 


0.167 


3.67x 10" 


-5 


0.104 



The uncertainty in the mass loss rate and dust-to-gas 
ratio derived from the dust shell is about 8% and 6% due 
to uncertainties in the assumed mass (varied to 0.6 M©) 
and vo (varied to km s _1 ). In the ideal situation the 
mass loss rates and dust-to-gas ratios from the two meth- 
ods would be identical to indicate a perfectly consistent 
model. The best overall agreement based on this analysis 
is achieved for a luminosity of 15 000 L©. 

Given all these considerations it is likely that the lu- 
minosity of IRC +10 216 is between 10 000 and 15 000 L© 
and its distance is between 110 and 135 pc. The present- 
day mass loss rate is (1.5 ± 0.3) xlO -5 M© yr _1 and the 
gas-to-dust ratio is about 700 ± 100. The dust opacity 
at 60 /im is found to be of order 250 cm 2 gr _1 . The CO 
abundance is 1.1 x 10~ 3 relative to H2. 

Compared to the previously proposed models men- 
tioned in the introduction this represents a somewhat 
lower mass loss rate, partly because these other models 
adopted larger distances as well in many cases. The CO 
abundance is larger than in most previous models, but in 
agreement with theoretical expectations. 
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The gas-to-dust ratio we find is larger that the canon- 
ical value of 200 often quoted. On the other hand there 
have been other recent studies that found larger gas-to- 
dust ratios for carbon stars as well, e.g. a gas-to-dust ratio 
of 550 for IRC +10 216 for a distance of 135 pc (Ivezic & 
Elitzur 1996), a value of 1230 for IRC +10 216 (Winters 
et al. 1997b), a value of 530 for AFGL 3068 (Winters et al. 
1997a). Note that in the model of Crosas & Mentcn (1997) 
the gas-to-dust ratio is fixed at 100, and they mention that 
a larger value would result in lines that are consistently 
too low. 
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